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Abstract. Current-voltage characteristics and dynamic conductance of the superconductor - 
constriction - superconductor junctions in GdO(F)FeAs polycrystalline samples with critical 
temperatures T^f?"''^ = 46 ^ 53 K were investigated. Two superconducting gaps, the large 
Al = 10.5 ±2 meV, and the small one As = 2.3±0.4meV were clearly observed at T = 4.2 K. 
The 2AL/fcBT,l;°'^='' = 5.5 ± 1 ratio gives support to the strong coupling mechanism which is 
responsible for the high Tc value. Temperature dependence of the large gap Al(T) indicates 
the presence of intrinsic proximity effect (in fc— space) between two superconducting condensates. 



Novel Fe-based superconductors [T] are now in the focus of intensive experimental research. 
Some of their features such as layered crystal structure and spatial separation of superconducting 
layers are similar to those of cuprates and MgB2. The stoichiometric 1 Ill-family compounds 
were shown to be antiferromagnetic metals with spin density wave (SDW) ground state |21I3]- 
Partial oxygen deficiency or fluorine substitution for oxygen induces superconductivity in the 
Fe-As layers. Rare-earth elements replacement also affects the critical temperature Tq. In 
particular, Tc for Gd-1111 oxypnictides can be lowered by replacement of Y for Gd [3] or 
increased by introducing Th instead of Gd The Tc = 56 K reported in Gdo.8Tho.2FeAs 
compound is today the highest one for Fe-based superconductors. 

Band-structure calculations showed the total density of states A^(0) at the Fermi level to 
be formed mainly by Fe Sd-states [H [71 [8]. A large Fe-isotope effect [9] and the correlation 
between Tq and A^(0) for different iron-based superconductors [10] give support to the phonon- 
mediated coupling in these compounds. The theoretically calculated Fermi surface for 1111- 
system compounds jTTl [121 [13 [H] comprises cylinder-like hole sheets centered at the T point 
and quasi-two-dimensional (2D) electron sheets at the M point of the reduced Brillouin zone. 

Here we present the superconducting gap measurements in GdO(F)FeAs by Andreev 
spectroscopy of ScS-contacts (superconductor - constriction - superconductor) using a "break- 
junction" technique. 

Nearly optimally doped Gd-1111 samples were prepared by high pressure synthesis described 
in detail in |151 I16j . Superconducting properties of the samples were tested by measuring 
temperature dependences of the AC-magnetic susceptibility and resistance R{T). Both showed 



a sharp superconducting transition at Tq = (52 53) K (the bulk value of the critical 
temperature was defined at the maximum of dR{T) / dT-cuive) . The two sets of polycrystalline 
GdO(F)FeAs samples were studied: "EF-series (GdOo.88Fo.i2FeAs, T^"''' = 53 ± IK) and 
"Kh"-series (GdOo.giFo.ogFeAs, T^^"^ = 52 it IK). Figure 1 shows typical temperature 
dependences of resistance for Kh-series samples. 

In order to determine superconducting gap values for GdO(F)FeAs, we used (i) Andreev 
spectroscopy [T7] of individual Sharvin-type [T5] ScS-contacts and (ii) intrinsic Andreev 
spectroscopy (intrinsic multiple Andreev reflections effect (IMARE) |19j . which usually exists 
due to the presence of steps and terraces at clean cryogenic clefts). ScS- Andreev nanocontacts in 
polycrystalline samples were prepared by the "break-junction" technique [2^ I21|. The samples 
of the typical size 1.5 x 0.5 x 0.1 mm^ were attached to the flexible sample holder with two 
current and two potential leads by In-Ga alloy that is liquid at room temperature. The holder 
with the sample was cooled down to 4.2 K, then a microcrack in the sample was generated by 
precise bending the sample holder. 

The resulting nano-contact is a mechanical contact of two clean cryogenically cleaved surfaces 
separated by a constriction; this geometry excludes impurity presence at the cryogenic clefts. 
The I{V)-, (i/(y)/(il^-characteristics of the contact were measured using a standard AG- 
modulation technique [22l [23] . 

According to the theory by Kiimmel et al. [23] for the SnS (superconductor - normal metal 
- superconductor) Andreev-type contact, the main features of its current-voltage characteristics 
(CVG) are the excess current at low bias voltages and the subharmonic gap structure (SGS). 
The latter shows series of dynamic conductance dips at bias voltages 
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due to the multiple Andreev reflections effect. For a multigap superconductor, several 
independent SGSs corresponding to each gap should be observed. The GVCs studied in this 
work are typical for the clean classical SnS-contact (with a constriction acting as a normal 
metal) with excess-current at low bias voltage, therefore, we believe, the theoretical model by 
Kiimmel et al. [21] is applicable to our break-junctions. Strictly speaking, the sharpest SGS 
are intrinsic of Andreev contacts of the highest quality with a small diameter a which is less 
than the quasiparticles mean free path I (i.e., in the ballistic limit) [18]. In such a case one may 
observe a large number of gap peculiarities, thus facilitating the gap energy measurements. In 
the case a ~ /, only few SGS minima would contribute to the dynamic conductance spectra, 
which in most cases corresponds to the studied GdO(F)FeAs ScS-contacts. 

In the framework of the "break-junction" technique and due to the layered structure of 
GdO(F)FeAs, the so-called stacks of Andreev contacts may be formed on cryogenic clefts; this 
type of junctions, the S — c— — S — c — 5-array, enables one to observe IMARE similar to 
that observed earlier in Bi-2201 [25] . For such an array, bias voltage for any peculiarities caused 
by the bulk properties at (i/(y)/dy-characteristic scales with number of sequentially connected 
contacts as compared with the single junction case. On the contrary, this is not the case 
for such CVG features that originate from any defects influence on the cleaved surfaces. One 
can determine by collating dynamic conductance spectra of Andreev arrays with different 
number of junctions in a stack and normalizing them to a single contact characteristics. The 
array contacts are a factor of N less sensitive to surface defects (which otherwise would seriously 
affect the properties of superconductor |26j). Thus, by making GVC measurements on the stack 
contacts, one can determine the true bulk gap values much more accurately. 

Figure 2 shows dynamic conductance of two ScS-Andreev contacts: single junction E13dl, 
lower black line, and an array Kh3d4, upper gray line. After normalizing the dI{V)/dV- 
characteristic of the array (A'^ = 2 for this case) to that for the single junction, we achieve a good 
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Figure 1. Normalized temperature de- 
pendence of resistance for polycrystalline 
GdOo.giFo.ogFeAs samples Khl (solid yellow 
circles), Kh2 (red line) and Kli3 (open blue 
circles) measured prior a microcrack forma- 
tion. The bulk critical temperatures Tq""^ = 
52 ± IK. 
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Figure 2. dl /dV-curves for the single- 
junction ScS-contact Elld6 (lower yellow line) 
and 2-j unctions Kh3d4 array (normalized to a 
single junction; upper red line). T = 4.2 K. 
SGS dips defining Al ~ llmeV and 
As ~ 2.5 meV are marked with tt-l and ns 
labels, respectively. 



coincidence of the spectra peculiarities. Such a procedure was done for all ScS-arrays, with N as 
a fitting parameter. The normalized dynamic conductance shows a set of dips at bias voltages 
V ~ 22 mV and 11 mV (marked by dashed vertical lines). The minimum at V ^ 7.3 mV is 
seen only at the Kh3d4-curve (marked by arrows) for stack contact. These three peculiarities 
may be associated with the SGS minima riL = 1, 2, 3, which leads to the estimated local value of 
the large superconducting gap Al ~ 11 me V. The minima at ~ it 5 mV at the ElldG-curve 
(marked by solid vertical lines) don't correspond to the expected bias voltage V ~ 7.3 mV for 
the 3rd harmonic = 3 of the large gap, and, hence, indicates the presence of the small gap 
As ~ 2.5 meV. The satellite minima above the til = 1 dips may originate from excitation of 
Leggett collective modes and require additional studies. 

For some contacts (such as E13 sample, single ScS-Andreev junction #dl), one can see the 
SGS originating from the small gap, but not only a single peculiarity. Figure 3 shows such 
series comprising ns = 1,2,3 at y ~ (4.4, 2.2, 1.3) mV (vertical dashed lines). Using eq. 
(1), we directly find the small superconducting gap value As ~ 2.2 meV. The additional fine 
structure seen in Fig. 3 at y ~ 5 mV may be caused by the As order parameter anisotropy and 
is subject of further studies. As the small gap SGS is located near zero bias voltages, it becomes 
smeared by the dramatic increasing of the dynamic conductance at the dl/dV-curve. Hence, 
the background subtraction or taking the second derivative of the CVC helps one to distinguish 
the small gap peculiarities, as was done to the spectrum at the Fig. 3. 

By a precise mechanical readjusting of the contact, we were able to observe characteristics 
of several ScS break-junctions within the same low-temperature experiment; the independent 
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Figure 3. I{V)- (blue line), dl/dV- 
(red line) and d^I/dy ^-characteristics (yellow 
line) for single ScS-contact E13dl. T = 4.2 K. 
The SGS dips associated with small gap 
Ag ~ 2.2 meV are marked by dashed vertical 
lines. Linear background is subtracted. 
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Figure 4. Normalized to a single junction 
SGS bias voltages VaL, Vns versus 1/n for 
several ScS-contacts (Tq = 46 53 K) 
investigated. The solid lines indicate the 
averaged gap values Al = 10.5 it 2 meV 
and As = 2.3 ± 0.4 meV at T = 4.2 K. 



local gap measurements improve the statistics of the results. To get values of Al,s and their 
scatter, in accordance with eq. (1), we plotted SGS minima positions Va as a function of their 
inverse number 1/n; this line should start at the (0,0) point. Our experimental data measured 
at T = 4.2 K are summarized on Figure 4 that shows the results for three samples with a 
similar Tq"^^ (Kh3d4 contact - down triangles, Elld6 - squares, E13dl - rhombs, E13d2 - circles, 
and E13dl3 - up triangles). All the data agree with each other and confirm the presence of the 
two distinct superconducting gaps. The averaged gap values obtained at T = 4.2 K \16\ I27j 
are: 

Al = 10.5 ± 2meV, As = 2.3 ± 0.4 meV 

for the critical temperatures T(?'^^^ = (46 -i- 53) K. The close agreement of the data for various 
samples shows the reproducibility of the measured gaps magnitude. 

IMARE spectroscopy technique also facilitates measurements of temperature dependence 
of the superconducting gaps. The subharmonic gap structure of the dynamic conductance of 
symmetrical SnS-contacts remains rather well-pronounced up to Tq. Importantly, to define the 
gap value no additional fitting is required [21]: at any temperatures up to Tq the gap A can be 
obtained straightforward using the expression (1). We note, that in case of alternative technique 
of the Andreev spectroscopy of SN-junctions, and for a two-gap superconductor, one needs to 
fit the resulting conductance curves using 7 fitting parameters. 

Figure 5 shows the dynamic conductance of ScS-array (sample Kh3, contact #d4, see Fig. 2) 
measured in the range 4.2 K < T < Tq. The curves are shifted vertically for clarity. The 
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Figure 5. The temperature affecting on the 
SGS (black labels) at normalized dynamic 
conductance of 2 ScS-junctions array Kh3d4 
(4.2 K < T < Tc). Al(4.2 K) ^ 11 meV. 
Local critical temperature Tq"^^^ ~ 46 K. 
The curves were shifted along the vertical 
scale for the sake of clarity. 
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Figure 6. Temperature dependence of the 
large gap AL(r) (circles) plotted under the 
data of Figure 5. r^°'=''^ ^ 46 Single- 
gap BCS-like behavior (dash-dotted line) is 
presented for comparison. Error bars show 
typical temperature smearing for Cooper 
pairs at two representative temperatures. 



peculiarities are washed out at Tq ~ 46 K. The "local" critical temperature is referred to 
as the temperature of the transition to the normal state in the contact area (with diameter 
a ~ Z). It is worthy to note that the Tq'^^^ can differ from the T^^"*, the averaged characteristic 
measured by non-local methods. The coincidence of the local Tq'^'^^ and the averaged Tq""^ may 
be observed only in the ideal clean and homogeneous sample. By taking the Al ~ 11 meV value 
(see Fig. 2) and Tq'^'^^ ~ 46K from Figure 5, for the Kh3d4 array we obtain an estimate for 
the BCS ratio, 2Al/A;bTc"^^^ f» 5.5; the ratio would however be reduced to 4.8 if we use the 
averag ed T^^"^ = 53 ± 1 K value, as in Refs. I27j . The former value 5.5 is close to that 
for LaOo.gFo.iFeAs (a lower- Tc analogue of Gd-1111) [28]. At the same time, the measured 
2AL/A;B7c°'^'^^-ratio exceeds the standard single-gap BCS value 3.52 for the weak coupling limit. 
Presumably, the 2Al/A;bTc value suggests the strong coupling in the condensates with a large 
gap value. As for the small gap, our data leads to 2As/A;bTc = 1.1 -i- 1.3 which is less than 
the BCS value 3.52. 

Figure 6 shows the temperature dependence of the large gap Al(T) obtained from our IMARE 
measurements (see Fig. 5) within the range 4.2 K < T < Tq"^^^. Similar measurements of the 
temperature dependence for the small gap represent a harder task, possibly require higher quality 
of samples, and have not been accomplished yet. The Al(T) dependence lies noticeably below 
the standard single-gap BCS-like curve (dash-dotted line in Figure 6). Note that the temperature 
smearing (see error bars plotted for the last data point) allows to resolve the gap value up to 
~44.5 K. The deviation in the Al(T) dependence could not be attributed to surface gap influence 



(i.e., it is not due to the real space proximity effect), because it is obtained from Andreev 
reflection measurements in the array contact, where the surface effects are reduced. Moreover, 
the AL(r = 0) energy coincides with the value for the single ScS-contacts (see Fig. 2,4). Were 
the real-space proximity effect be responsible for the deviation, the Al(T = 0) would differ; our 
experiment, however, shows this is not the case. Such behavior resembles the Ao-(T)-dependence 
for MgBa [SHI EQl EI] and Al(T) for FeSe l23]. Theoretical studies [321 [331 El] explain such a 
deviation by interband coupling effect. By analogy, we assume that similar deviation for Gd- 
1111 may arise from nonzero interband coupling. The reduced BCS-ratio for the small gap 
also supports this conclusion and suggests that the "weak" superconductivity in the low gap 
condensate may be induced by fc-space (internal) proximity effect between two superconducting 
condensates ^5]. According to Refs. [521 E31 El], the observed shape of the A(r) curve (see 
Figure 6) is typical for the "driving" condensate in the presence of the second (driven) condensate 
with a small gap. 

The existence of the large superconducting gap (2Al/A;bTc""^ > 4) in the 1111-family 
compounds with critical temperatures 40 K < Tq < 53 K was confirmed by tunneling 
spectroscopy ("break-junction" technique) [361 EZ], point-contact Andreev reflection (PCAR) 
spectroscopy of SN-junctions [381 ESI HOI HH [121 [131 [H], and angle-resolved photoemission 
spectroscopy (ARPES) [l5]- To the best of our knowledge, there are no other data available for 
Gd-1111 to compare with. Thus, comparing our results for GdO(F)FeAs with the aforementioned 
data on other 1 Ill-superconductors, we find rather good agreement with 2Al/A;b2c values 
determined by "break-junction" technique \36\ I37j and in some PCAR-measurements |43t I44j. 

In conclusion, we studied the dynamic conductance of ScS-contacts in GdOi_a;Fa;FeAs 
{x = 0.09, 0.12) polycrystalline samples (Tq"^'^^ = 46 53 K) by ScS Andreev spectroscopy, 
in the temperature range 4.2 K < T < Tq- We detected two superconducting gaps, the 
large Al = 10.5 it 2meV, and the small Ag = 2.3 ± 0.4 meV. The observed temperature 
dependence of the large gap AL(r) suggests the existence of fc-space proximity effect between 
two superconducting condensates and, therefore, nonzero interband interaction. The estimated 
2Al/A:bTc"^'^^ = 5.5ibl ratio exceeds the BCS limit 3.52 and indicates the strong coupling regime 
in the "driving" condensate. 
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